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It is noted that if large iV c continuum reduction applies to an observable, then that observable is 
independent of temperature for all temperatures below some critical value. This fact, plus the fact 
that mesons and glueballs are weakly interacting at large N c is used as the basis for a derivation 
of large N c continuum reduction for the chiral condensate. The structure of this derivation is quite 
general and can be extended to a wide class of observables. 
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One of the most intriguing recent developments in 
QCD is the prospect of large N c continuum reductionflj, 
Hi El El Underlying this is the two decades old idea of 
Eguchi and Kawai (EK) that as N c , the number of de- 
grees of freedom in QCD, goes to infinity, the relevant 
spatial size in the dynamics goes to zero 5]. This was 
originally formulated on the lattice and it was shown that 
the Wilson loop was computable on a single plaquette, 
from which it is deduced that the free energy density 
is computable on a single plaquette. Unfortunately, al- 
though the EK analysis is formally correct, it has been 
of limited practical significance there is a phase transi- 
tion separating the EK reduced theory from the physical 
point as one approaches the continuum limit 0. How- 
ever, recently it was suggested that while the continuum 
theory cannot be reduced to a point — it can be reduced 
to a Euclidean box of finite size (in terms of physical 
units) while keeping the exact large N c infinite volume 
result 0,0, 0,0]. The guts of this argument is that pro- 
vided the box size is larger than the critical size for which 
the phase transition occurs, the EK reduction works and 
one can simultaneously avoid the phase transition and 
thus connect with the continuum limit. Numerical sim- 
ulations in three and four dimensions^, 0, provide 
strong, but not compelling, evidence that this scenario is 
correct. The phenomena is referred to as large N c con- 
tinuum reduction. 

It has also been conjectured that the same phenomena 
of large N c continuum reduction occurs for the chiral con- 
densate: namely, the chiral condensate as computed in a 
finite box at large N c becomes independent of box size 
for boxes larger than some critical valueQ. The conjec- 
ture is motivated by consideration of random matrix the- 
ory (RMT) which is believed to become universal for the 
lowest modes of the Dirac operator; this presumed uni- 
versal behavior is in terms of the natural scaled variable 
in RMT, Zfc = Afc-/V C L 4 E where a is the chiral condensate 
and Afe is the fcth eigenvalue of the Dirac operator. One 
expects the same universal behavior by taking L — > oo 
or N c — ► oo. 

This conjecture may seem paradoxical: it is well known 
that spontaneous symmetry breaking cannot occur for 
systems of spatial extent. The resolution to this paradox 
is simple: a necessary (though not sufficient) condition 



for spontaneous symmetry breaking is not infinite spa- 
tial size but rather an infinite number degrees of freedom 
(which, of course, occurs for a system of infinite spatial 
extent). There can be an infinite number of degrees of 
freedom in a system of finite spatial extent if one has 
an infinite number of internal degrees of freedom; this 
is precisely what happens for QCD as N c — ► oo. There 
is numerical evidence that the chiral condensate can be 
computed on a finite box by an extrapolation to infinite 
N c ; moreover, the evidence is consistent with the possi- 
bility that the value of the condensate is independent of 
the box size beyond some critical size 7J. If the contin- 
uum reduction holds, then one can do reliable large N c 
simulations on relatively size boxes leading to relatively 
inexpensive computations. This prospect is exciting; reli- 
able and converged lattice QCD results may be available 
in the near term — albeit for a colorful world. 

The present letter has two main purposes: to point out 
that the continuum reduction has very strong implica- 
tions for the phenomenology of large N c QCD at nonzero 
temperature and to provide a basis in well-known large 
N c phenomenology for understanding why continuum re- 
duction must occur for the chiral condensate (and for the 
expectation values of a wide class of other operators). 
These two issues are intimately related. 

First consider the phenomenological implications for 
large N c QCD at T ^ 0. The central result is that any 
quantity which is subject to the continuum reduction at 
large N c has the property that for any temperature less 
than some critical temperature T c , the quantity is un- 
changed from its T — value. This result is not exact at 
finite N c but any deviations from the T = are sublead- 
ing: 
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where Q(T) is some quantity subject to large N c con- 
tinuum reduction. Thus, for example, if the chiral con- 
densate is subject to the continuum reduction (as con- 
jectured in ref. 0) then in a large N c world as the tem- 
perature of the system as increased from zero, the chiral 
condensate will remain unchanged from its T = value 
until the temperature reaches a critical point. 
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The physical basis is quite simple. The logic of large 
N c continuum reduction implies that when the Euclidean 
space box size is larger than the critical size, the results 
at large N c do not depend on the size or shape of the 
box: if all of the sides of the box are larger than the 
critical length (L c ), then the quantity does not depend 
further on the lengths of the sides. Thus for a box with 
three of the sides of infinite extent and the fourth of fi- 
nite size (greater then L c ), then a quantity subject to 
continuum reduction does not depend on the length of 
this fourth side. However, Euclidean field theory with 
one side of finite length and standard boundary condi- 
tions (periodic for bosons, anti-periodic for fermions) is 
just finite temperature field theory with T = 1/L (where 
L is the length of this side). Thus, any quantity for which 
large N c continuum reduction holds will be independent 
of temperature for T < T c where T c — l/L c up to 1/N C 
corrections. 

This behavior is consistent with what is known about 
the temperature dependence of the chiral condensate. 
Consider the chiral condensate at fixed low temperature 
in the combined large N c and chiral limits. For sim- 
plicity, consider the case where the chiral limit is taken 
prior to the large N c limit; this avoids ambiguity about 
defining the chiral condensate by limiting the sensitivity 
to the ultraviolet. Moreover, with this ordering chiral 
perturbation theory at finite temperature can be used to 
compute an analytic expression for the variation of the 
chiral condensate with T (at zero quark mass) . Up to 
order T 6 the result is: 
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where (qq)T is the chiral condensate at temperature T, 
is the pion decay constant (evaluated in the chiral 
limit) and A q is a parameter which depends on higher- 
order counter terms in the chiral lagrangian. Recall the 
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standard N c scaling of j-n'-J-n ~ N c ■ Using this scaling 
and taking temperatures of order iV°, and combining it 
with eq. (0 yields the result that 
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There is no variation of the chiral condensate with tem- 
perature in a large N c world up to the order calculated, 
as is required by large N c continuum reduction. 

Equation y| is highly suggestive and its important to 
understand whether it holds generally. Fortunately, it is 
possible to show on very general physical grounds that 
the chiral condensate is independent of temperature at 
sufficiently low temperatures in a large N c world. Ob- 
viously, this is of phenomenological interest (at least in 
a large N c world). In addition, this observation can be 
used as the physical basis of a demonstration that the 
phenomena of continuum large N c reduction applies to 
the chiral condensate. Moreover, the logic underlying the 
derivation can be generalized to the expectation values 
of a wide class Lorentz scalar local operator and not just 
to the this chiral condensate. 

The key to this general derivation is the fact that stan- 
dard 't Hooft counting|9( implies that at large N c , mesons 
and glueballs have masses of order iV° and are weakly in- 
teracting. V^ ncs0 n — meson ^ N c , V^ ncson _gl uc ball ^ N c 

and V^i U ebaii-giuebaii ~ N~ where V indicates the 
strength of the two-body interaction (four-point func- 
tion). Now consider an extended system in which the 
density of mesons and glueballs is of order N®. In such 
a system the thermodynamics will be simply that of a 
gas non-interacting mesons and glueballs plus corrections 
which vanish at large N c . Thus, at temperature of order 
iV°, the density of a given species, s, of a meson or glue- 
ball is given by the standard expression for noninteract- 
ing bosons: 



/d k 1 
— 3 p(T, m s , k) (1 + OiN- 1 )) with p(T, m s ,k) = 7==- (4) 

where d s = (2J S + 1) (2/ s + 1) is the degeneracy factor for a species with angular momentum J s and isospin / s . Note 
in terms of N c counting p s (t) ~ AT?. The free energy density as measured relative to its zero temperature value is 
computed in a similar manner and is given by 

G r = -Tj2d s J ^log(l- e -^^^/ ^ )(l + 0(iV- 1 )) (5) 
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where the summation is over species of boson. The sub- convention one sets Q (0) to zero and this distinction is 
script r indicates that the free energy is measured relative unimportant. However, in the present context we will 
to its T = value: Q r {T) = Q{T) — £7(0). Normally, by differentiate with respect to an external source, such as 
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the quark mass, and G(0) has nonvanishing dependence 
on these sources. 

Before proceeding it is useful to recall that strength 
of the interaction between a meson or a glueball and a 
baryon is of order unity in N c counting and thus is not 
weakly interacting. So, if the system has a baryon density 
of order unity, eqs. (@J and will be spoiled. However, 
the baryon mass, unlike the meson or glueball masses, is 
order N c 0] and hence, for temperature of order AT°, the 
density of baryons will be of order e~ N " and is negligible 
in the large N c limit. 

Consider the chiral condensate, i.e., the expectation 
value of qq where both quarks are understood to be at 
the same space-time position. Of course, this operator is 
not well defined without specifying some prescription for 
how to regulate it at large momentum. But, as we will sec 
for the present purpose, the details of this prescription 
turns out to be largely irrelevant. For the moment let us 
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neglect this issue and proceed formally. To simplify issues 
associated with renormalization it is useful to discuss the 
RG invariant operator mqq (where m is the quark mass) 
rather than qq itself. From the standard form of the 
grand partition function as a Euclidean space functional 
integral it is immediately apparent that 



m(qq) 



T dlog(Z) 
V dm 



(6) 



where the system is confined to a three-space volume V 
(taken to infinity in the thermodynamic limit) and is of 
length = 1/T in the time direction; m is the quark 
mass. Defining the free energy density in the usual way 
as Q = —y\og(Z) it is apparent that (qq) — J^. Com- 
bining this with eqs. iQJ andJSJ along with the definition 
of Q r yields: 
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where <r s the sigma commutator for the species s and is 
defined by a s = m^^-. 

The key point for the purposes here is the N c scaling. 
By conventional large N c scaling arguments 0, E3 it is 
straightforward to see that 
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AT" 1 for s e {glueballs} 
iV c ° for T ~ 7V C ° 
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Combining these scaling rules with eq. J7|) yields, 



m(qq) T - m(qq) T =o 
m(qq) T =o 
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This is precisely of the form of eq. Q as required by large 
N c continuum reduction. 

A few words about renormalization are in order. In 
practice the only viable way to numerically estimate the 
functional integral for the partition function is via the 
lattice and one might worry that using a lattice regu- 
larized theory might invalidate eq. ©. However, if one 
uses fermions which respect chiral symmetry ( eg. overlap 
fermions 0) then one can use eq. © to define a regu- 
lated chiral condensate in a consistent manner. As noted 
above the combination mqq appears in the Lagrangian 
and as such has no operator mixing and no anomalous 
dimension. However, this does not mean the expecta- 
tion value is defined unambiguously. Indeed, even for 



a free field theory this quantity diverges in the ultravio- 
let. To render the operator finite it is typically defined at 
some scale p where fluctuations at scales above p are sub- 
tracted off. Provided p is large enough this is computable 
perturbatively. Moreover in the regime A >> p >> T, 
where A is the cutoff scale of the regularized theory, (eg. 
1/a for a lattice regularized theory) then eq. 10 can be 
used to compute m(qq)T — m(qq)T=o and the result is 
independent of the choice of p. The reason for this is 
simple — the same high momentum fluctuations are sub- 
tracted from both the finite T and T = cases so that 
their difference does not depend on the subtraction point. 
Thus, the numerator of eq. © is independent of p and 
of order N®. The denominator does depend on p but for 
any value of p, it is of order A^. This yields the scaling 
in eq. © for any p. Finally it should be noted that in the 
chiral limit, the ratio in eq. © is independent of p; both 
the numerator and denominator have /i-independent con- 
tributions at leading chiral order m with all p dependence 
coming in at order m 2 or higher. 

It is very easy to understand the origin of eq. 
Equation has a very simple interpretation. It rep- 
resents the amount each hadron contributes the spatially 
integrated operator mqq weighted by the density of the 
hadrons. Clearly the contribution from a single meson is 
order N® as the operator is a one-body quark operator 
and the number of quarks in a meson is independent of 
Af c . Since the density of mesons is also of order unity 
one sees that the shift in the condensate from T = is of 
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order unity. In contrast the condensate scales with the 
number of colors since it effectively counts the number 
of active quarks species. The scaling of eq. © follows 
from this. From this qualitative discussion, it should be 
clear that although the preceding derivation was for the 
chiral condensate, the same basic feature should occur 
for generic local Lorentz scalar operators. The expecta- 
tion value at T — will be altered at finite T due to 
the contributions from non-interacting mesons and glue- 
balls and such contributions are characteristically smaller 
in N c counting. Of course, the issues of operator renor- 
malization for more complicated operators becomes more 
involved. 

The temperature independence is only supposed to 
hold below the phase transition; if the particles don't 
interact one might worry that no phase transition can 
occur. However, as N c — ► oo QCD has an infinite num- 
ber of narrow mesons and glueball and the number of 
accessible states below some energy E can grow expo- 
nentially with E. If the number of accessible meson and 
glueball states below E goes as e~ R / Tc , then the total 
density of mesons and glueballs ^ s p s (T) will diverge as 
T — > T c . This is just the Hagedorn phenomenon and 
its description predates QCD[l3; presumably, the QCD 
phase is increasing well described as a Hagedorn transi- 
tion as N c — > oo. Of course, as the Hagedorn transition 
is approached, the 1/N C suppression in meson-meson in- 
teractions is overwhelmed by an enhancement due to the 
diverging density; one can no longer neglect interactions 
and the derivation presented here breaks down. Thus at 
or above T c , nothing prevents the chiral condensate (or 
other observable) from differing from its T = value at 
large N c . 

The phenomenological demonstration that below T c 
the chiral condensate (and other observables) are inde- 
pendent of T at large N c can be used to derive the full 
result of continuum reduction. Note, that virtually noth- 
ing in the derivation depended on the spatial size being 
infinite. If one repeated the derivation for a the system 
quantized in a box with sides of length L x ,L y and L z 
and standard boundary conditions (periodic for bosons, 
anti-periodic for fermions) , the only alterations would be 
in eqs. 10} and © for which the following substitutions 
would be made: 




The key point is that with this substitution eq. J7J) con- 
tinues to hold as do the scaling rules in eq. (|SJ| and thus 
so does eq. ©■ 

Now imagine one starts with a box of infinite extent 
in all four directions. Labeling one of these directions 
as time one can use the arguments leading to eq. © to 



reduce the length of this dimension without altering the 
chiral condensate (to leading order in 1/N C ) providing 
we do not pass through the phase transition; this is just 
a temperature shift. Next one can relabel the axes and 
denote another axis as the time direction and reduce it. 
This also simply becomes a temperature shift on a box 
which is now finite in one spatial dimension. However, as 
was just argued, temperature independence of the chiral 
condensate holds for the case of boxes with finite sides. 
One can similarly reduce the other two sides leading to 
the result that at large N c the chiral condensate for an 
arbitrary size finite box is equal to that of an infinite box 
provided no phase transition is crossed. This is large N c 
continuum reduction. 

The logic underlying this derivation of large N c contin- 
uum reduction for the chiral condensate is complimentary 
to that of ref. 7]. It is based only on well known phe- 
nomenology of hadronic interactions and simple features 
of finite-temperature field theory and thus avoids a re- 
liance on assumptions about the applicability of random 
matrix theory for this context. However, the present ar- 
gument is quite limited compared to that in ref. 0] in at 
least one critical way. In contrast to ref. it provides 
no insight into why chiral symmetry is broken in the first 
place; it merely explains why the chiral condensate does 
not depend on box size (at large N c ) for boxes larger than 
a critical size. 
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